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Abstract 

We find, and analyse, the exact solution of two friendly directed walks, modelling 
polymers, which interact with a wall via contact interactions. We specifically consider 
two walks that begin and end together so as to imitate a polygon. We examine a 
general model in which a separate interaction parameter is assigned to configurations 
where both polymers touch the wall simultaneously, and investigate the effect this 
parameter has on the integrability of the problem. We find an exact solution of the 
generating function of the model, and provide a full analysis of the phase diagram 
that admits three phases with one first-order and two second-order transition lines 
between these phases. We argue that one physically realisable model would see two 
phase transitions as the temperature is lowered. 



1 Introduction 

The adsorption of polymers on a sticky wall, and confined between two walls, has been the 
subject of continued interest pTJ El El SI El El El El El QUI QUI 02] - This has been in part due to 
the advent of experimental techniques able to micro-manipulate single polymers [13} [bfl [15] 
and the connection to modelling DNA denaturation p^|IT7lfT8l[T9|[20l[2T|l22l[23]. 

Consider a polymer in a dilute solution of good solvent, so that it is in a swollen state 
|21]. If such a polymer is then attached to a wall at one end the rest of the polymer 



drifts away due to entropic repulsion. On the other hand, if the wall has an attractive 
contact potential, so that it becomes 'sticky' to the monomers, the polymer can be made 
to stay close to the wall by a sufficiently strong potential, or for low enough temperatures. 
The phase transition between these two states is the adsorption transition. The high 
temperature state is desorbed while the low temperature state is adsorbed. This pure 
adsorption transition has been well studied JTJ [2j [25j El I2S] exactly and numerically, and 
has been demonstrated to be second-order. 

There has been recent interest [H] in ring polymers, modelled by self-avoiding polygons, 
being adsorbed onto the walls of a two-dimensional slit. In that work, models in which both 
sides of the polygon could interact with each of the walls were considered. This provides us 
with one of the motivations for the model here, where we consider two directed walks that 
begin and end together, so forming a polygon. We consider directed walks because they 
often admit exact solutions, while the more realistic self-avoiding walks do not. Moreover, 
we consider such a pair of walks interacting with a sticky wall, allowing different interactions 
when one or both walks are near the wall. To allow for a simple realisation of the model we 
consider so-called friendly directed walks (rather than the ubiquitous vicious walks) where 
the two walks may share edges of the lattice. However, we do not allow the walks to cross 
and so there is always an upper walk/polymer and a lower walk/polymer. 

Other physical motivations for two-walk models have appeared in the literature. In 
particular, one may model DNA-denaturation in this way — for example the Poland- 
Scheraga models [271 EE] • It would be interesting to see if the techniques described below 
could be used to find exact solutions of the DNA unzipping transition in the presence 
of an adsorbing wall, such as that discussed in [29]. To do so we would add a contact 
interaction between the two walks, rather than the double- visit interaction discussed here. 
A manuscript on this topic is in preparation [30J . 

As we investigated the model another motivation for its interest became apparent; the 
full two parameter model is not amenable to one of the standard methods of solving multiple 
walk models. The Lindstrom-Gessel-Viennot lemma [31, E2] (which was also considered 
earlier in a probabilistic context by Karlin and McGregor [33J) decomposes the solution of 
models of multiple vicious walks into combinations of single walk problems. The lemma 
implies that the generating function of a multiple walk model would be governed by a D- 
finite (Differentiably-Finite) function, but the solution of our model is not D-finite. Despite 
this, we are still able to solve the model. 

We have solved our model in two ways. Firstly, we use the obstinate kernel method (see 
[31] for an overview of the technique) to give a formal solution of a functional equation as 
a constant term formula. This constant term can then be evaluated explicitly. Secondly, 
we also use a 'primitive piece' decomposition that allows us to give an explicit solution in 
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terms of hypergeometric-type sums. 

Our solution allows us to fully analyse the model and we find a rich phase diagram. 
In particular, we find three phases that meet at a special point. There are three phase 
boundaries; two are second-order and one is first-order. Intriguingly, we find that arguably 
the most physically realisable one-parameter case of our model would have two phase 
transitions on lowering the temperature. 

In the next section (Section]!]), we formally define our model. In Section[3j we formulate 
a functional equation obeyed by an extended generating function and provide a 'constant 
term' solution in Section [4] via the obstinate kernel method. We provide an alternate 
explicit solution in Section [5] that illustrates an underlying structure in the solution which 
arises combinatorially. In Section |6j we analyse the phase structure and phase transitions of 
the model. In the final section (Section [7]) we discuss the functional nature of the solution 
and summarise our results by recasting them in terms of some physical parameters of a 
family of single parameter models. 



2 Model 

We consider a pair of directed walks above a wall on the upper half-plane of the square 
lattice, taking steps (1,1) or (1,-1). These walks may touch but not cross; such walks are 
sometimes called friendly walks. Further, we consider those pairs of walks that begin at 
the point (0, 0) and have equal length. Let <p be a pair of such walks and the set of all such 
walks be Q. We define \(p\ to be the length of the walks. 

To these configurations we add an energy — e a for visits of the bottom walk only (single 
visits) to the wall, and an energy — when both walks visit a site on the wall simultaneously 
(double visits), excluding the first vertex of the walks. The number of single visits to the 
wall will be denoted m a ((p), while the number of double visits will be denoted rrid(<p). 

Later in the paper we will specialise to those configurations, <p in which both walks 
start and end on the wall. Since every such configuration has at least one double visit (the 
final vertex), we have m^ip) > 1. The partition function for our model is 

Z n (a,d)= e ima ^ y£a+md ^ y£d)/kBT (2.1) 

(p^ \(p\=n 

where T is the temperature and ks the Boltzmann constant, and associated Boltzmann 
weights are denoted a = e £a//fesT and d = e £d ^ kBT . The thermodynamic reduced free energy 
of our model is given in the usual fashion as 

n(a,d) = lim n -1 log (Z n (a, d)) . (2.2) 
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A configuration of length 10 in our model with single and double visits marked appears in 
Figure [T] 




Figure 1: Two directed walks of length 10 of our model that begin and end on the surface. 
There are two single and two double visits marked. The left-most (start) vertex of the two 
walks on the wall is not counted clS db double visit. 

To find the free energy we will instead solve for the generating function 

oo 

G(a, d;z) = J2 Z n(a, d)z n . (2.3) 

n=0 

The radius of convergence of the generating function z c (a, d) is directly related to the free 
energy via 

k(q, d) = log(z c (a, cfT 1 ). (2.4) 

3 Functional Equations 

To find G, we consider walks (p in the larger set, where each walk can end at any possible 
height. Let us first consider a — d — 1. In this case we construct the expanded generating 
function 

F(r,s;z) =F(r,s) = ^z^Ms^ 2 , (3.1) 

^ eft 

where z is conjugate to the length \ip\ of the walk, r is conjugate to the distance [ip\ of the 
bottom walk from the wall and s is conjugate to half the distance \tp~\ between the final 
vertices of the two walks. Since the distance between the endpoints of the walks changes 
by or ±2 with each step, and the endpoints start together, it is always an even number. 
Further, we let [r- 7 s fc ]F(r, s) denote the coefficient of r J s fc in the generating function F(r, s). 
We use [rijFfrjs) to denote the coefficient of r 7 in F(r,s) which is a function of s and 
similarly [s k ]F(r, s) gives a function of r. 
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Let us now return to general a and d. All pairs of walks can then be built using a 
standard column-by-column construction. Translating this into its action on the generating 
function gives the following functional equation 

F(r, s) =1 + z (r + - + - + -) • F(r, s) 
\ r r s J 

-z(- + -Y \r°]F(r, s)-z-- [s°]F(r, s) 
\r r J s 



+ z{a - 1)(1 + s) ■ [r x ]F(r, s) + z(d — a) ■ [r s ]F(r, s). 



(3.2) 
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Figure 2: Adding steps to the walks when the walks are away from the wall. There are 
four possibilities. 

We explain each of the terms in this equation. 

• The trivial pair of walks of length gives the initial 1 in the functional equation. 

• Every pair of walks may be extended by appending directed steps to their endpoints 
in four different ways (see Figure |2]) . 



Top walk 


Bottom walk 


Generating Function 


(1,1) 


(1,1) 


r • F(r, s) 


(1,1) 


(1,-1) 


- r -F(r,s) 


(1,-1) 


(1,1) 


L s -F(r,s) 


(1,-1) 


(1,-1) 


\-F{r, s) 



Appending steps in this way may result in the bottom walk stepping below the wall 
(Figure [3]). Thus, when the bottom walk is at the wall, we cannot append any steps 
that will decrease the height of the bottom walk. These forbidden configurations are 
counted by 



1 s 
: I - + - 



• \r°}F(r,s) 



: 1 I + J ) ' F{0.s). 



(3.3) 
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Figure 3: The first boundary term in the functional equation removes the contribution 
from the walks that are produced by appending a SE step to the bottom walk when its 
endpoint is on the wall. 




Figure 4: The second boundary term in the functional equation removes the contribution 
of walks that cross. Such configurations are produced when one appends steps to walks 
that end at the same vertex as shown. 



• Similarly, at no time can the top walk pass below the bottom walk (Figure [4]). Thus, 
if the two walks are touching, we forbid the distance between them to decrease. These 
configurations are counted by 

z- ■ [s°]F(r,s) = z- -F(r,0). (3.4) 

• This accounts for the possible pairs of walks without the interaction parameters. We 
can now incorporate the interaction parameters. In order to do this, we have to add 
in all walks we want to mark with a and subtract the non-weighted version of those 
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exact same walks from the model (Figure [5] — left). In order for the bottom walk to 
touch the wall, it must be at height 1 initially and then step down (with no restriction 
on the top walk). Hence we get the term 

z(a-l)(l + s) ■ [r^Ffas). (3.5) 




Figure 5: Configurations that lead to single (left) and double (right) interaction terms. 



A similar method can be used to incorporate d into the model (Figure [5] — right). 
One step before both walks touch the wall they will both be at height 1. All such 
walks have already been accounted for when incorporating a into the model and so 
must be replaced. This results in 



z(d-a) ■ [rV].F(r, s). 



(3.6) 



The functions [r 1 ]F(r, s) and [r 1 s°]-F(r, s) can be simplified in terms of F(0, 0) and F(0, s). 
By extracting the coefficient of r°s° in the functional equation, we obtain 



F(0,0) = l + zd- [rVji^s). 



(3.7) 



At a combinatorial level, this states that a pair of walks that end at the wall is either a 
trivial configuration or obtained by appending a pair of SE steps to the end of a pair of 
walks that end at height 1. Similarly, we can extract the coefficient of r° in the functional 
equation to obtain 



F(0, s) = 1 + za(l + s) ■ [r x ]F(r, s) + z(d — a) ■ [r s ]F(r, s). 



(3.1 



This has a similar combinatorial interpretation to the previous case. These equations can 
then be combined to simplify the functional equation 



1 



1 s r 
r+ -+-+ - 
r r s 



■ F{r, s) 



+ 1 



1 

a 



zs 
r 



F(0,s) 



zr 



F{r,0) 



F(0,0). (3.9) 
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We will use the above form of the equation in what follows. The polynomial coefficient on 
the left hand side is called the kernel K(r, s; z) = K(r, s), and its symmetries play a key 
role in the solution: 



K(r,s) 



I S T 

l-z[r + - + - + - 
r r s 



(3.10) 



4 Solution of the functional equations 

In what follows we use the obstinate kernel method. The discussion below is self-contained, 
but we refer the reader to the paper of Bousquet-Melou and Mishna [31] for a general 
description of this technique. 

4.1 Solution of the functional equations when a = 1 



When a = 1, the functional equation (3.9) simplifies to 



K(r, s) ■ rsdF(r, s) = rs - zsd (1 + s) ■ F(0, s) - zr 2 dF(r, 0) + rs (d - 1) F(0, 0). (4.1) 

We use the kernel method which exploits the symmetries of the kernel to remove boundary 
terms (ie the functions F(r, 0) and F(0, s)) in the above equation. The kernel is symmetric 
under the following two transformations: 

r \ /s 



(r,s)^[r,-), (r,s)^-,sj. (4.2) 

These transformations generate a family of 8 symmetries (sometimes referred to as the 
'group of the walk' — see |34j) 



M, I r, r -), C-, 4) , IV (I IV (l 4\ ( T -_, r -\, and f t S ) . (4.3) 



s y Vr r z / \s s J \r s J \r r z J \s s J \r 

We make use of 4 of these transformations — those which only involve positive powers 
of r. To be precise, 

K(r, s) ■ rsdF(r, s) = rs - zsd (1 + s) ■ F(0, s) - zr 2 d ■ F(r, 0) + rs (d - 1) • F(0, 0); (4.4a) 



r 



Kir,— ) F r, — = 1 + — -F 0,— - zr^d • F(r, 0) + — (d - 1) • F(0, 0); 

s/sys/s s \ si \ s I s 



(4.4b) 

, r r" \ dr" I r r" \ r" zar' I r 2 \ ( r 2 \ zr 2 d ,r r 3 , 

R (?Tj-^(?Tj = ^-^( 1 + Tj- F ( 'Tj-^- F ^ Q ) + ^^- 1 )- i? ( ' ) ; 

(4.4c) 

(r 1\ d- r p r F (q Z ^ r2 F ( T o) + ^ _ ^ F(0 0) 

\s's/s 2 \s' s J s 2 s \ s) \ S J s 2 \s' J s 2 

(4.4d) 
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All of these transformations were chosen so that the kernel remains unchanged, and so that 
the substitution only involves positive powers of r. We can then eliminate the boundary 
terms by taking an alternating sum of the above equations: 



Eqn(4.4a) - Eqn(4.4b) + Eqn(4.4c) - Eqn(4.4d). 



(In the case where a 7^ 1, a similar method holds, except that then we must multiply some of 
the equations by non-trivial coefficients to eliminate boundary terms.) After simplification 
we obtain 

K (r, s) ■ (linear combination of F) = 

r(s — l)(s 2 + s + 1 — 1 



s 



2 



-^(l + (d-l)F(0,0)) 

- zd(l + s)sF(0, s) + Zd{l + S) F (0, . (4.5) 

We can now remove the left-hand side of the equation by choosing a value of r that sets the 
kernel to zero — provided all the F's on the left-hand side remain convergent. The kernel 
has two roots and we choose the one which gives a positive term power series expansion in 
z with Laurent polynomial coefficients in s: 



r is: z ) = r 



l _ /1 _ 4 ( 1±£ )^ 



2(1 + s)z 

K ' n>0 



h -L e ^2n+l„2n+l 

+ . («) 



where C n = ^;( 2 ") is a Catalan number. This is chosen so that K(f, s) = 0, and so that 
all the various substitutions are convergent. More precisely, since f = O(z), the functions 
F(f, s), F(r, f 2 / s), F(r/s, r 2 /s) and F(r/s, 1/s) are all formally convergent power series in 
z with Laurent polynomial coefficients in s. 

We are not able to use the other root of the kernel (with respect to r) since it is 0(z~ l ). 
If we were to substitute this into the functional equation, then F(r, s), F(r, r 2 / s), F(r/s, r 2 / s) 
and F(r/s, 1/s) would not converge within the ring of formal power series. This follows 
since the coefficient of z n in F(r, s) has degree n in r and so substituting r 1— > 0{z^ 1 ) will 
map terms in this polynomial to all powers of z including the constant term. 

When we make the substitution r^fwe can rewrite the coefficients of the right-hand 
side so as to not explicitly involve z — since now z — (f + 1/f + f/s + s/r)~ l . 

r(s — l)(s 2 + s + 1 — f 2 ) , . , , „. NN dfs 2 . 1 dr „/ 1 

= J )A 5 1 1 + (d - 1 )F 0, 0)) - — - F 0, a) + F 0, - 

(4.7) 

Because we are primarily interested in F(0, 0) — the generating function of pairs of 
walks that start and end on the wall — it is convenient to rewrite the equation so that there 
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are no powers of s or f in the denominator of the coefficients and so that the coefficients 
of F(0, s) and F(0, 1/s) are independent of f. 

1' 



= cfe 4 F(0, s) - dsF (^0, -j-(s-l)(s 2 + s + l- f 2 )(s + f 2 ) (1 + (d - 1)F(0, 0)) . 

(4.8) 

Consider the coefficient of s 1 in the above equation, or rather by dividing the equation 
by s consider the constant term, that is the coefficient of s° in the equation. This leads us 
to calculate F(0, 0) effectively as a constant term in the variable s. For ease of calculation 
and display we will continue with calculating the coefficient of s 1 in equation (4.8). 

Since -F(0, s) is a power series in z with polynomial coefficients in s, the term ds 4 F(0, s) 
does not contain any coefficients of s 1 . Similarly, F (0, -) is a power series in z with 
polynomial coefficients in s~ x , so the term dsF (0, i) contributes only dF(0,0). For the 
remaining term, we consider the coefficient [s 1 ](s — l)(s 2 + s + 1 — f 2 )(s + f 2 ). Expanding 
the expression and then collecting the exponents of f gives: 

(s - l)f 4 + (1 - s + s 2 - s 3 )f 2 + s(l - s)(s 2 + s + 1). (4.9) 

We need to consider the expansion of f 2 and f 4 . Lagrange inversion [35] gives: 

r(s; z) 



°° r* h -I- e^n+l 



n=0 



/•(■s:,) 2 = X: ^ ^- iob ) 

n=0 



f(s; z) 4 



E _4_/2n + 4\(l + .)^ w 

^2n + 4V n / s n K J 



n=0 



and, more generally, 



n=0 v 7 

Computing the coefficient of a particular power of s in f 2 or f 4 reduces to finding the 
coefficient of powers of s in (1 + s) n s~ k which are just binomial coefficients: 

[a i]( a -l)^ = f;-^J>^^; (4.11a) 

n=0 

ro-^- ; -A^E ( „ 6 g ( ;; 3) ^; 

n=0 7 

[s 1 ]s(l-s)(s 2 + s + l) = 1; (4.11c) 
[s%s -l)(s 2 + s + l-r 2 )(s + r 2 ) = l + jr, ^^ + ^ Cg^ a . (4.11d) 
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Hence extracting the coefficient of s 1 in equation (4.8) gives 

= - ( 1 + £ . yS/ r) • (1 + (d - 1)F(0, 0)) - d • F(0, 0). (4.12) 
Solving the above when d — 1 gives 

g(i.M) = i+E ( „S ( t^ ^. <" 3 > 

and hence for general d we have 

F(0,0) = G(M ;,) = d + (1 G ^ M , z) . (4.14) 

In Section [5] we will see that the algebraic structure of this solution that gives G(l,d;z) 
in terms of G(l,l;z) arises naturally from a combinatorial construction. Moreover, this 
structure extends to the a ^ 1 case. 

4.2 Solution of the functional equation when 1 

The general a, d case can be solved by the method applied above, however, it is sufficient 
to study the case d = a which can be resolved more cleanly. As mentioned above the 
algebraic structure that allows G(l, d; z) to be expressed in terms of G(l, 1; z) extends to 
give G(a, d; z) in terms of G(a, a; z). We shall see that explicitly in Section 15} When d = a 



the functional equation (3.9) simplifies to 

K(r, s) ■ F(r, s)a 2 rs = (ar — r — za — zas)as ■ F(0, s) — zr 2 a 2 ■ F(r, 0) + ars. (4.15) 

The symmetries we used above can be reused to remove boundary terms. As above we 
take an alternating sum of transformed equations, but now we must multiply some of the 
equations by a non-trivial factor chosen to eliminate boundary terms. The left-hand side 
becomes 

LHS = a^K(r, s) (sF(r, s) - (r, + g, - *F g, i)) , (4.16) 

where 

zas — ars + rs + zar 2 , „ _ „. 

L = — . (4.17) 

zas — ar + r + zar z 

The right-hand side simplifies to 
RHS = as 2 (l + s- a)F(0, s) + a(l + s - as)F (0, l/s) 



(r 2 + s)a(s - l)(ar 2 + as - 2s - s 2 
ar 2 — r 2 — s 



(4.18) 
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Again, we have attempted to massage the functional equation into a form in which the 
coefficients of F(0, s) and F(0, 1/s) are independent of r. Unfortunately, we cannot com- 
pletely clear the denominator of the above functional equation, and we found it simplest 
to work with the above expression. 

Following the method used in the a = 1 case, we can eliminate the left-hand side further 
by choosing a value of r that sets the kernel to 0. We choose the root which gives a positive 
term power series expansion in z with Laurent polynomial coefficients in s. Recall that f 
is given by 

(1 + s ) 2n+1 z 2n+1 



£C. " T ",„ ' . (4.19) 



r{s; z) =r — 

n>0 

Substituting r t- > f eliminates the left-hand side of the functional equation and we again 
consider the coefficient of s 1 in the resulting right-hand side. Again, this can be converted 
into a constant term expression for our generating function. The term as 2 (l + s — a)F(0, s) 
does not contribute to s 1 . The term a(l + s — as)F (0, contributes a(l — a)F (0, 0). 
For the remaining term, we consider the expansion of the expression as a series in a. The 
coefficient of a 1 is 

[a ^(r 2 + s)(s-l)(af 2 + as-2s-s 2 -l) + 



ar 2 - — r 2 



and hence 



[a i s , f(r 2 + s)(s-l)(ar 2 + as-2s-s 2 -l) = _^ 
ar 2 — r 2 — s 



Higher powers of a are (for k > 1) 



Il a(f 2 + s)(s — l)(af 2 + as — 2s — s 2 — 1) 



ar 2 — f 2 — s 



(s — l)(f — l)(f + l)(f — s)(f + s) / f 



r.2 



f 2 \s + f 2 

(S " 1)S2 -(s + l)%s-l)).(^-) k (4.22) 



\s + 1) 2 Z 2 1 / \S + f 2 y 

fc 



To extract the coefficient of s 1 , we need to consider the expansion of J in z. This 
exponential term simplifies, and we can use equation (4.10d) to obtain 



r 



2 \ k /1 1 „\ fc 



k ,i + A \ k _^- *> (2 P + a (g + i) 2 ^+ 2fc 2fe 



= z * f fc = > — — ' " 1 ~; — (4.23) 
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We will expand the two terms in (4.22) individually. For the first term, we get 



[s + l) 2 z 2 \s + f 2 

k fk + 2p\ f f2k -2 + 2p\ /2fc - 2 + 2p 



E 

p>0 



fc + 2p V p 



k + p — 1 



k + p-2 



z 



2fc-2+2p 



(4.24) 



We can extract the coefficient of z 2n from the above equation by making the substitution 
n = k — 1 + p. We obtain 



s + 1) 2 Z 2 \ s + f 



;2 \ fc 



fc 



2n - fc + 2 \ n + 1 



2n - k + 2 



2n 
n 



2n 
n — 1 



(4.25) 



Therefore 



(s + 1) 2 Z 2 \s + f 



s2 \ k 



E 

n>fc-l 



fc 



2n - fc + 2 V n + 1 



2n - fc + 2 



2n 
n 



2n 
n — 1 



„2n 



(4.26) 



Following a similar argument for the second term in equation (4.22), we have 

-2 \ k 



[<?](, + l)°(s - 1) ■ 

= E 



p>0 



s + f 2 

fc fk + 2p\ f f2k + 2p + 2\ f2k + 2p + 2 
k + p 



k + 2p\ p 



k + p — 1 



Making the substitution n = k + p, we get 



z 2 V](s + l) 2 (s-l) 



r 



2 \ k 



S + f z 



fc (2n — k 



2n — k\ n 



2n + 2\ (2n + 2 



n + 1 



We can then substitute the summation over p with a summation over n. 



[s^s + l) 2 (s - 1) 



r 



2 \ fc 



S + 7" 



E 

n>k 



k (2n — fc 



2n — k\ n 



2n + 2 

71 + 1 



n 



2n + 2 



. (4.28) 



.2» 



(4.29) 



When n = fc — 1 in the above equation, the summand reduces to when fc > 2. So it is 
possible to adjust the range of the summation by adjusting for the fc = 1, 2 cases separately. 
In those cases, the combined correction terms are a 2 and — 4a 3 z 2 respectively. Thus, we 
can rewrite 

52 \ k 



[£\{8 + \)*{8-\) 



r 



s + f z 



E 

n>k-l 



2n — k V n 



2n-k 



2n + 2 
n + I 



2n + 2 

n 



z 2n , (4.30) 
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with known correction terms for k = 1,2. Combining these summands, we get that for 
n > k — 1: 



k 



2n-k + 2\ n + 1 



2n-k + 2 



2n 
n 



2n 
n — 1 



k (2n — k 



2n — k\ n 



2n + 2 
n+1 



2n + 2 
n 



k(k + 1) (2 + An - kn - 2k) (2n - k\ (2n 



(k-l-n)(n + l) 2 (k-2n)(n + 2)\ n 



n 



(4.31) 



Thus taking the coefficient of s 1 when r = f in equation (4.18) and accounting for the 
correction terms, we get 



= a(a - 1)F(0, 0) - a + a 2 - Az 2 a 3 

k(k + 1)(2 + An- kn - 2k) (2n-k 
(k-l-n)(n + l) 2 (-2n + k)(n + 2) \ n 



E aM E 



k>\ n>k 

We exchange the order of summation to give 



2n 



n 



z ln . (4.32) 



2„3 



= a(a - 1)F(0, 0) - a + a 2 - Az 2 a 



n+1 



E 2 ' 2 "E 



k(k + 1) (2 + An - kn - 2k) (2n - k\ (2n 



n>0 



'(& — ! — n){n + l) 2 (-2n + k)(n + 2) 



n 



n 



,fc+i 



(4.33) 



The extraction of the coefficient [a k z 2n }F(0, 0) requires rearranging the a(a — 1) coeffi- 
cient in front of F(0, 0). We express the above equation as: 



n+1 



= a(a - 1)F(0, 0) - £ z 2n £ Q n , fe ,a fc ' +1 



(4.34) 



n>0 fc'=l 

for some integers (5 n ,fe'- This can be rearranged to give: 

n+1 



i71 (°»°) = -(E ;z2 "E^" / 

v n>0 fe'=l 

n+1 

£ 2 ™ Qn,k>" 
v n>0 fe'=l 



A:' 



k' 



1 — a 

E' 

vfc">0 



.fe" 



(4.35) 
(4.36) 



The coefficient of a fc from the above is summation of all contributions from k' and k" such 
that k' + k" = k. Thus: 



n+1 



F(0,0) = ^^^a fe ^g„„ fc , 



(4.37) 



n>0 fc=l fc'=l 
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In other words, extracting the coefficient [a z 2n ]F(0, 0) requires a summation of a finite 



number of the Q n ,k' terms which is obtained from (4.33). 

fc 



r„* J»i F( n n) = V k'(k' + 1)(2 + An ~k'n- 2k') (2n - k'\ (2n 
^ Q (k' -l-n)(n + l) 2 (-2n + k')(n + 2) \ n 

k(k + l)(k + 2) (2n-k\(2n 
~ (2n-k){n + l) 2 {n + 2) \ 



n 



n 



We finally obtain 

F(0,0) = G(a,a) 



V- 2n^ a k k(k + l)(k + 2) (2n\(2n-k 
^ {n + mn + 2){2n-k)\n) \ 



n>0 k=0 



n 



(4.38) 
(4.39) 



(4.40) 



This agrees with results due to Brak et al. [361 EZ] for a closely related model obtained 
using very different method — we discuss this more fully in the section [5} 



5 Alternate solution 

An alternate technique for finding the generating function relies on factoring the pairs of 
walks at each double- visit. First, let us define 

G(a, d) = G(a, d; z) = F(0, 0; a, d; z). (5.1) 

We will frequently hide the z dependence for convenience. Breaking up our configurations 
into pieces between double visits gives 

G(a, d; z) = - — !— (5.2) 

1 — dPya] z) 

where P(a; z) is the generating function of so-called primitive factors. This is quite analo- 
gous to the classical factorisation of a single Dyck path. These primitive factors are pairs 
of friendly Dyck-paths which contain no double-visits to the surface other than their first 
and last vertices. Rearranging this expression gives 

p , x = G ( a > d;z)-l = G(a, a; z) - 1 
dG{a, d; z) aG(a, a; z) 

This last expression allows us to calculate P(a; z) from a known expression for G(a, a; z) 
- such as that given in the previous section. Alternatively, one could use results from 
previous work by Brak et al. [361 137] . In those works, the authors considered a vesicle 
model which corresponds exactly to the case d = a — their vicious walk model can be 
transformed into the friendly walk model considered here, by moving the upper vesicle 
boundary down by 2 units. 
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Brak et al. use the Lindstrom-Gessel-Viennot lemma [3TJ [32] to express the partition 
function of the pair of walks in terms of the partition function of a single walk. Namely, 

r mnf \ S 2n {\)S 2n+i (a) - S 2n+2 (l)S 2n+2 (a) 

[z \&{a, a) = 5 , (5.4) 

cr 

where S 2n (a) is the partition function of a single Dyck path of length 2n above a wall, and 
a is conjugate to the number of visits, ie 

S 2n ( a ) = [z 2n ] 2(2-a + aVl- 4z 2 ) (5.5) 



n + fc + lV^ - ^/ 2n — A; \ w 

fc=0 1 1 \ / fc = i v / 

where this last formula is taken from [36]. When a = 1 we recover the well-known Catalan 
number result, and a well-known central binomial result when a = 2: 

5 2n (l) = C n = -4r ( 2n> ) , S 2n (2) = f 2n ) . (5.7) 

In light of these simple expressions one can write G(a, a; z) as double sum of products of 



binomials. Using equation (5.3) we write G(a, d; z) in terms of G(a, a; z) 



n( a\ aG(a,a) 

G(a,d) = - — -, ——. r, (5.8) 

v ' d+(a- d)G{a,a) K ' 

where 

oo 

G(a, a) = a' 2 ^ [C n S 2n+4 (a) - C n+1 S 2n+2 (a)] z 2n , (5.9) 

ra=0 



which simplifies to the expression in equation (4.40) found in the previous section. 



5.1 Solutions at a = 0, 1 and 2 

Since the partition function S 2n (a) takes simple values at a = 0, 1, 2, we have 



G(l, 1; z) — [C n C n+2 — C 2 +1 ] 

n=0 

= V 12(2n + 1) ( 2n Xz 2 - (5 10) 

oo 

G(2,2;z) = ^C n C n+1 z 2 ™, (5.11) 



n=0 
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and 



lim 

a-s>0 



G(a, a; z) 



E 



n=l 



12(2ri - 1) 



n 2 (n+ l) 2 (n + 2) V n-1 



2n - 2 



= z*G(l,l;z). 



(5.12) 



We can use these together with equation (5.8) to derive expressions for G(l,d) (that 

G(2,d] 

Eoo 
n=0 



agrees with (4.13) and (4.14)) and G(2,d) by simple substitutions. That is 

12(2n+l) 



G(l,d) 



(2n\ 2 2n 

(n+l) 2 (n+2) 2 (n+3) \ n) 41 



d+(l-d)EZ 



12(2n+l) 



(n+l)2( n+ 2)2( n +3) \ n 



(5.13) 



and 



G(2,d) 



2i( 



A little further work also gives 

1 1 

G(0,d) 



l-dz*G(l,l) l-dz^Z -„n# 



12(2n+l) 



/"2rA 2 2n ' 



(5.14) 



(5.15) 



(n+l) 2 (n+2) 2 (n+3) 

This last expression can be derived combinatorially by noting that in the limit a — > 
single visits are forbidden. In this limit, the primitive pieces are in bijection with all walks 
counted by G(l, 1); any primitive piece can be transformed into a pair of walks counted by 
G(l, 1) by moving them 1 lattice unit up and gluing edges at the start and end. 



6 Analysis of phase structure and transitions 
6.1 Phases 

We now turn to the phase diagram of the model which is dictated by the radius of con- 
vergence of G(a, d; z) as a power series in z. Denote the radius of convergence by z c (a, d). 



Equation (5.2) shows that the singularities of G(a, d; z) are those of P(a; z) and the simple 



pole at 1 — dP(a;z) = 0. Denote this latter singularity by Zd(a,d). Equation (5.3) shows 
that the singularities of P(a; z) are related to those of G(a, a; z) which are known from 
[361 ETJ. 

In particular, the radius of convergence of G(a, a; z) is 

a < 2, 

p{a) = r.o.c G(a, a; z) = ^ r (6.1) 

z a (a) a > 2. 




2d 



For a < 2, the thermodynamic phase is related to zj, and is the desorbed phase in which 
the walks drift away from the wall and the mean number of visits is 0(1). When a > 2, 
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z a (a) dominates and the lower walk adsorbs onto the wall and the number of visits is 0(n). 
At a = 2, there is a second-order phase transition and a jump discontinuity in the specific 
heat (the second derivative of the free energy). In both of these phases, the upper walk 
drifts away from the wall, and the number of doubly- visited vertices is 0(1). 

In the full a-d model there are 3 phases, two of which are described in the previous 
paragraph. In the third phase, associated with the simple pole at Zd(a,d), we shall see 
that the number of doubly- visited vertices is 0(n). In what follows, we name these three 
phases associated with Zb, z a and zj, desorbed, a-rich and d-rich respectively. 

6.2 Desorbed to a-rich transition 

In [52], it was shown that the asymptotic behaviour of the singular part of G(a,a; z) near 
its radius of convergence is given by 

4z) 4 log(l -Az) a<2, 

4z) 2 \og(l-Az) a = 2, (6.2) 

1 /2 

■ z/z a (a)) a > 2, 

where z a (a) = v/ °~ 1 . It is important to notice that for all a, the singularities are convergent 
and therefore G(a, a; z) is convergent on its radius of convergence p(a). 

If we fix d at some small value, and then increase z from towards p(a), then P(a; z) 
increases from to P(a; p(a)). Since d is small, and P(a; p(a)) is finite, 1 — dP(a; p(a) > 
and so the only singularities of G(a, d; z) will be those of P(a; z) and so those of G(a, a; z). 

Thus for small values of d there is a phase transition on moving a through 2 which 
describes the transition from the desorbed phase to an a-rich phase as occurs in [36J. This 
adsorption transition has been well-studied previously and is unusual in that it has a jump 
discontinuity in the second derivative of the free-energy rather than a divergence. 

6.3 Desorbed to d-rich transition 

Let us restrict our attention to a < 2 and consider the effect of increasing d. The argument 
in the previous subsection breaks down as soon as 1 — dP(a; p(a)) = 0. Call this value 
d c (a) = P(a; p(a))~ 1 . 

Fix a < 2 and d > d c and consider increasing z from towards p(a) = 1/4. The 
function P(a; z) is an increasing function of z (since it is a positive term power series) and 
so it increases towards P(a; p{a)). Since d > d c , P(a; z) will reach the value d~ l before it 
reaches P(a,p(a)) and the simple pole will occur when z = Z4, where Zd is the solution of 

P(a; Zd(a, d)) = gT 1 (6.3) 



G(a, a; z) 



A_(l- 
A (l- 
A4.fl 
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and Zd < Zb = 1/4 in this region. 

Hence, for a < 2, there is a phase transition where Zd and ^ coincide at 

, , s „/ . . i aG(a, a: 1/4) 

4 a = Pa; 1/4 - 1 = r , 1 ''I V 6.4 

G (a, a; 1/4) — 1 

In order to determine the density of the singly- and doubly-visited vertices in the d > d c 
phase consider the partial derivatives of Zd(a, d) with respect to a and d. Since Zd(a, d) 
is defined by dP(a; z d (a,d)) = 1, the derivatives of Zd(a,d) with respect to a and d are 
non-zero and so there are positive densities of both singly- and doubly-visited vertices. 

Now let us turn to the order of this transition; this can be determined by examining the 
behaviour of P(a;z) close to z = 1/4 which is determined by the behaviour of G(a,a;z) 



see equation (5.3). Close to z = 1/4 we can write 

G(a, a; z) = G ana i y ti c (o,', z) + G S i ngu i ar (a; z), (6-5) 



where the behaviour of G S i ngu i ar is given by equation (6.2). Consider an expansion of 

G analytic about 2 = 1/4 

G analytic^ z) w G(a, a; 1/4) + d(l - 4z) H (6.6) 

for some non-zero constant C\. The linear correction dominates the dominant singular term 



in G S i n g U i aT . Expanding equation (6.3) about z = 1/4 gives 



i- - = P(a; 1/4) -P(a;z d ), 
d c d 

d -^^ Pl (l-Az d ) (6.7) 
d 2 c 

for some non-zero constant p\. Hence there is a linear relationship between the location of 
the d-rich singularity, Zd, and the distance from the phase boundary. Since the free-energy 
of the system is — log Zd, this also implies the free-energy in the d-rich phase changes linearly 
with d. On the other hand in the desorbed phase where z^ dominates, the free-energy is a 
constant. From this we see that there is a jump discontinuity in the first derivative of the 
free-energy and hence this is a first-order transition. Note that the above argument will 
also work mutatis mutandis at a = 2. 

We can observe at finite length a characteristic bimodal probability distribution in the 
number of doubly- visited vertices — see Figure [6] 

6.4 a-rich to d-rich transition 

The analysis of the previous section can be adapted to the case a > 2 with some important 
differences. The transition is driven by the singularities z a (a) = associated with 
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Figure 6: A plot of the probability that a conformation of length 128 has k doubly- visited 
sites at a — 1, d — 10.3. This value of d corresponds to the approximate location of the 
peak in the specific heat at this length. 

single-visit adsorption, and the singularity Zd(a, d) associated with double-visit adsorption. 



Again, Zd(a,d) is the solution of equation (6.3). These two singularities coincide when 
d = d c (a) given by 



d c (a)=P[a;^- = , V . 7 ■ (6.i 



2a J Gr^^.y^T] _ , 



Turning to the order of this transition, we again decompose G(a, a; z) into its analytic 



and singular parts. Observe that close to z a (a), G S i ngu i ar , given by equation (6.2 ), dominates 
the linear part of G ana iy tic . Hence we deduce that 

d c (a) - d ^ p 2 f — — z d (a, d) J , (6.9) 

z d (a, d) as 1 +p 3 (d- d c (a)) 2 (6.10) 

2a 

for some nonzero constants P2,P3- Therefore the free-energy has a jump discontinuity in 
its second derivative on varying d across the transition, and this is a second-order phase 
transition. This is very similar to the desorbed to a-rich transition. 
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6.5 Phase diagram 



We have established that there are 3 thermodynamic phases; desorbed, a-rich and d-rich. 
We remind the reader that m a (<p) and md{(p) denote the number of single and double visits 
of (p. 

If we define 

(m a ) 



A(a, d) = lim 



and 



n— >oo n 

then in the desorbed phase we have 



V(a,d)= lim {md) 



n— >oo fl 



A = V = 0, 

while in the a-rich phase we have 

A > and 

and in the d-rich phase has both 

A > and 

In Figures M and M we plot and respectively. 



256 



v = o, 



V>0. 




(6.11; 



(6.12) 



(6.13) 



(6.14) 



Figure 7: A plot of the density of a visits calculated at length n = 256. This highlights the 
region where it tends to a non-zero constant and corresponds well to the regions where z a 
and Zd dominate. Note that for fixed a and increasing, large d we expect that the density 
of a visits decreases though remains positive. 

The phase boundary between the desorbed and a-rich phases occurs at 

a = 2 for d < 4(2). (6.15) 
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m d 
256 




Figure 8: A plot of the density of d visits calculated at length n = 256. This highlights the 
region where it is tends to a non-zero constant and corresponds well to the region where 
we have shown that Zd is the dominant singularity. 

Note that this phase boundary is, unsurprisingly, independent of d. We can compute d c {2) 



exactly using the results of Section 5.1 



64 

G(2,2; 1/4) = Y,CnC n+l lQ- n = 8 - — ; (6.16) 



3tt' 

n>0 



2G(2,2) 16(8 - 3tt) , 

* (2) = = " 1L55159579 - (6 ' 17 > 



In a similar way we can compute d c (0) and d c (l): 

30 

= i«e ~ 14.81234030; (6.18) 
1657T — 512 

= 8(512-165.) ^ ( lg) 

v ; 4096 - 1305tt v ; 

The transitions to the <i-rich phase from the desorbed and a-rich phases are given by 

, . . „. , -i aG(a,a; p(a)) 

d c (a) = P a; p(a))- 1 = — 1 P \ )} 6.20 

G(a, a; p(aJJ — 1 



where p(a) is given by equation (6.1). We plot P(a; p(a)) and G(a, a; p(a)) in Figure [9j 

In the limit as a — >• oo, d c (a) — > 2a + o(a). As a — > oo, the generating function G(a, a) is 
dominated by those configurations which have a maximal number of visits. In this case, the 
lower walk simply zig-zags along the wall and the upper walk is effectively unconstrained 
by the lower. Hence 

lim G(a, a; z) = V C n a n z 2n = 1 - ^ l ~ AaZ \ (6.21) 
a-^oo ^-^ 2az 2 

n>0 
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Substituting z = p(a) then gives G(a,a; p(a)) = 2 + o(l) from which the asymptotics of 
d c (a) follows: 

<i c (a)~2a as a — > oo. (6.22) 
In the next section we compute the above asymptotic form in more detail. 



2 I ii 0.1 








10 20 10 20 

a a 



Figure 9: A plot of both G(a,a; p{a)) and P(a; p(a)). The dotted line indicates a = 2. In 
the plot of G(a, a; p(a)) we have also marked the asymptotic form computed below. 

Combining all of this information gives the phase diagram for our model, which we 



present in Figure 10 It is interesting to note that the three transition lines meet with the 
two critical lines forming an angle. Classically this would indicate mean-field like behaviour 
of a bicritical point. If true, this mean-field behaviour would be interesting to understand. 



6.6 Asymptotics of the d-rich-a-rich phase boundary 

We now consider how the ci-rich a-rich phase boundary, d c (a), behaves for large a in more 
detail. From equation (5.4), we write [z 2n a k ]G(a, a; z) in closed form 



k(k + l)(k + 2) (2n\ (2n - k 



z a ]G(a, a: z) = -, —r-, — -^— / — - — — 

J v ' ' ; (n + l) 2 (n + 2)(2n- k) \ 



n \ n 



(6.23) 



We seek the asymptotic form of G(a, a; p(a)) as a — > oo, and so we need to evaluate the 
asymptotics of 

. ^ (a -If (2n\^k(k + l)(k + 2) (2n-k\ k 

G(a, a; p(a)) = g (w - ^ (Jg { Jl k) \ n f- ^ 

Expanding this slightly further gives 

. .. v-^ 1 [2n\ (-iyk(k + l)(k + 2) fn\f2n-k\ i+k 
Q; P ^ =^ (n + lY{n + 2)^ U ) g (2n - k) %){ n f 



n>0 v ' v ' v ' j,k=0 



(6.25) 
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Figure 10: The phase diagram of our model. The three phases are as indicated and the 
first-order transition is marked with a dashed line, while the two second-order transitions 
are marked with solid lines. The three boundaries meet at the point (a, d) = (2, 11.55 . . . ). 

Substitute j = £ — k to get 

= v - 1 (2n\ y> e __ 2n y, (-l) l - k k(k + l)(k + 2) ( n \f2n-k 

~ ^ (n + l) 2 (n + 2)4 n \n J ^ a ^ {2n - k) \£-k)\ n 



So now the coefficient of a is 

.o 1 = y-_ 1 

+ 1 V n 



n=0 



while 



la' 1 



la- 3 



n=2 



and 



la~ A ] 



(6.26) 



K 2 ] = (6.29) 
^ (n + l)(n + 2) Vn/ v ; 



f 3(n-2) (2n\ 5 

^( n +l)(n + 2)UJ 4 1 j 
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These simple forms continue as far as we have observed. This then gives 

G(a, a- P (a)) ~ 2 -\ + ^ + 1 + A + ^ + + 0(0 . (6 . 32) 

We can then plot this against our numerical estimates of G(a,a; p(a)). One should note 
that the series G(a, a, p(a)) converges very slowly for a > 2. Since we know the summands 
decay as n~ 3 ^ 2 , we can assume that the partial sums, s n , grow as A + Bn~ x l 2 . We can then 
accelerate the convergence of the series by estimating A with the sequence 

A n = s n (n + yfn\Jn - 1) - s n _i(n - 1 + y/ny/n - 1). (6.33) 

This combination was chosen by solving the pair of simultaneous equations s n = A + 
Bn- 1 ' 2 , s n _! = A + B(n - 1)~ 1/2 . We found that sequence A n converged far faster than 
the partial sums s n . 



7 Discussion 



7.1 Nature of solution 

In the section [6] we demonstrated that when a — 1, the model undergoes a phase transition 
at d = d c {l) = 4o 5 9 6 2 r^o^ ■ Since this is not an algebraic number, it follows that the 
generating function of the model does not satisfy a linear differential equation in z with 
integer polynomial coefficients in a, d and z. That is, it cannot be D-finite. 

Consider, to the contrary, that the generating function G(l, d; z) = [r°s°]f(r, s; 1, d; z) 
is a D-finite power series in z with integer polynomial coefficients in d. By definition, it 
satisfies a non-trivial linear differential equation of the form 

d k G dG 
Pk(d;z)-^ + ...p^z)— +p G(l,d;z) = 0, (7.1) 

where the Pj(d; z) are integer polynomials in d and z. By standard results in the the- 
ory of linear differential equations, the singularities of G(l, d; z) are zeros of the leading 
polynomial Pk(d; z). 

For small d we know that the dominant singularity of G(l, d) is Zb — 1/4. At the critical 
value of d, there is a change in dominant singularity from Zf, to Exactly at the critical 
value, Zf, = Zd = 1/4. Thus the discriminant of Pk{d; z) with respect to d must be zero 
at this point. Since the discriminant is a polynomial in z, d with integer coefficients and 
z = 1/4, it follows that this critical value of d must be an algebraic number. Above we 
showed that d c {l) = ^Hzl^l which is not algebraic and thus G(l, d; z) is not D-finite. 
A standard result [HE] on D-finite series states that specialisation of D-finite series are 
themselves D-finite and thus G(a, d; z) cannot be D-finite and nor is f(r, s; a, d; z). 
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The Lindstrom-Gessel-Viennot lemma combines the partition functions of single-walk 
models — equivalent to sums and Hadamard products of the underlying single-walk gen- 
erating function which is algebraic (this is true quite generally — see |39j). Any finite 
combination of Hadamard products and sums of algebraic or D-finite generating functions 
remains D-finite (38] and thus the Lindstrom-Gessel-Viennot lemma (alone) cannot be 
applied to decompose the model considered here into single-walk problems. 

That being said, the Lindstrom-Gessel-Viennot lemma can be combined with a factori- 
sation argument to yield a solution as we demonstrated in section |5j 

7.2 Fixed energy ratio models: r-models 

Finally, let us now consider the family of physical models parameterised by — oo < r < oo 
where 

Ed = re a and so d = a r (7.2) 

that allows us to summarise our results. Let us call these r-models. 

For any r-model the high temperature phase is the desorbed state. The model effectively 
already analysed by Brak et al. [361 E2j has d = a and so r = 1. In this case there is a 
single low temperature phase being the a-rich phase. Given there are no additional phase 
boundaries for a < d one can deduce that for all r < 1 the model goes from the desorbed 
state at high temperatures through a single second-order phase transition to the a-rich 
phase at low temperatures. 

The special point in our phase diagram where (a, d) = (2, 11.55. . .) where the three 
phases meet occurs in the r-model with 

_ log(11.55...) 
r = r t = — = 3.53.... (7.3) 

For all r > r t there is a single low temperature phase which is the d-rich phase: the 
transition on lowering the temperature is now first-order. 

Since we have shown above that d c (a) ~ 2a as a — > oo one can now argue that for all 
1 < r < r t the r-model has two phase transitions on lowering the temperature. At very 
low temperatures the model is in a d-rich phase while at high temperatures the model is in 
the desorbed state. At intermediate temperatures the system is in an a-rich phase. Both 
transitions, from desorbed to a-rich, and a-rich to d-rich, are second-order transitions with 



jump discontinuities in the specific heat. In Figure 11 we plot the fluctuations in a-visits 
as a function of temperature at length 128 for the r = 2 r-model: two peaks occur in these 
fluctuations. 
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Figure 11: A plot of the fluctuations in the number of a-visits, m a , for length n = 128 as 
function of a clearly showing two peaks. 

If one argues that a physically realisable model is one where both walks pick up the 
same energy when they touch the surface together then the model is the one with r = 2, 
that is d = a 2 . It is interesting to see that this model contains two phase transitions: one 



at a = 2 and the other at a « 3.301 found by solving equation (6.20) for d c = a c . In any 
case, we have a family of adsorption models that have one or two low temperature states 
and which the order of the transition changes as the parameter is varied. We have analysed 
this model using an exact solution and fully delineated its behaviour. It will be of interest 
to analyse the behaviour of this model in a slit. 
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